An efficient matrix method for coupled systems of variable fractional order differential equations

We establish a powerful numerical algorithm to compute solutions of coupled system variable fractional order differential equations. Our numer?ical procedure is based on Bernstein polynomials. The mentioned polynomials are non-orthogonal and have the ability produce good results as compared some other method like wavelet. By differentiation integration, operational matrices formed. On using obtained matrices, proposed reduced algebraic Using MATLAB, we solve given equation for required results. Graphical presentations maximum absolute errors illustrate Some useful features our sachem those that need no discretization or collocation technique prior develop matrices. Due these computational complexity much more reduced. Further, efficacy enhanced by increasing scale level. also compare with Haar wavelet justify useful?ness adopted method.